Non-perturbative approach to correlations in two-dimensional vortex liquids 



Joonhyun Yeo and M. A. Moore 
Department of Physics, University of Manchester, Manchester, MIS 9PL, United Kingdom. 

(February 1, 2008) 



in 

OS 



in 



> 
o> 



in 
o 
in 

o> 



i 

C 

o 
o 



We calculate the renormalized quartic vertex function of the Ginzburg-Landau model for a super- 
conducting film in a magnetic field by summing an infinite subset of diagrams, the so-called parquet 
graphs. Using this non-perturbative solution, we obtain the structure factor of the two-dimensional 
vortex liquid. We find growing crystalline order in the system as the temperature is lowered. Our 
results suggest that the length scale characterizing the crystalline order diverges only in the zero- 
temperature limit, which indicates the absence of a finite-temperature phase transition to the vortex 
lattice phase. 
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The effect of thermal fluctuations in type II supercon- 
ductors in the presence of a magnetic field has been a 
focus of theoretical and experimental interest since the 
discovery of high T c superconductors. Mean field theory 
predicts the existence of the mixed phase where the mag- 
netic field penetrates the sample in the form of flux lines 
which form a triangular array known as the Abrikosov 
vortex lattice [jjj. Thermal fluctuations are believed to 
be responsible for the melting of the Abrikosov vortex 
lattice into a vortex liquid state. The melting of the 
two-dimensional vortex lattice has been studied in asso- 
ciation with the dislocation unbinding mechanism lead- 
ing to a continuous transition However, the question 
whether the two-dimensional vortex liquid undergoes a 
phase transition at all into the low-temperature ordered 
state is still controversial. The high temperature pertur- 
bation expansions of the Ginzburg-Landau (GL) model 
within the lowest Landau level (LLL) approximation have 
been analyzed for evidence of the transition to the or- 
dered phase at low temperatures [gj. Some numerical 
simulations j| of the same model revealed a weak first- 
order transition between vortex liquid-solid phases. On 
the other hand, the Monte-Carlo (MC) simulation per- 
formed in Ref. ||ona spherical geometry showed no sign 
of a finite temperature transition. 

In this Letter, we present a non-perturbative analysis 
of the GL model for a superconducting film in a magnetic 
field within the LLL approximation. We make a resum- 
mation of an infinite number of diagrams summing all the 
parquet graphs ||. We find a clear indication of grow- 
ing crystalline order in the two-dimensional vortex liquid 
as the temperature is lowered. However, within our ap- 
proximation, we find no evidence for a finite-temperature 
phase transition. In fact, our results are consistent with a 
zero-temperature scaling argument Q according to which 
the length scale characterizing the spatial order in the 
vortex liquid diverges only in the zero temperature limit. 

A different kind of non-perturbative scheme for the 
same model has been employed by Tesanovic and co- 



workers to calculate various thermodynamic quanti- 
ties. This method, however, contains an adjustable pa- 
rameter which is the generalized Abrikosov ratio, (3 a- 
In the present nonperturbative analysis, where no ad- 
justable parameters are used, all the physical quantities 
including Pa are calculated as functions of temperature 
in a straightforward manner. 

The starting point of our analysis is the Ginzburg- 
Landau free energy for a superconducting film in a per- 
pendicular magnetic field B = V x A given by 
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where a, (3 and m are phenomenological parameters. For 
a thin film superconductor, fluctuations in the vector po- 
tential A can be ignored. We restrict the order parameter 
^ to the space spanned by the LLL wavefunctions. This 
approximation is valid near the upper critical field H C2 . 
In the symmetric gauge, where A = V> x )> the LLL 
is fully described by an arbitrary holomorphic function of 
the variable z = x + iy multiplied by an exponential fac- 
tor: 9(x,y) = exp(—^-z*z)(j)(z), where fx = y / e*B/~h is 
the inverse magnetic length. In the LLL approximation, 
the GL free energy becomes 



F[(j>] = / dz* dz{anc 



-l 2 |0(z)| 2 + ^ 2 W>(z)| 4 }, 

(1) 

where J dz*dz = J d 2 r and an = a+Tie* B /2m — is the 
point where the mean field instability occurs. The tem- 
perature will be represented by a dimensionless parame- 
ter, ctT = an ^/^n/ Pfi 2 . In the high temperature pertur- 
bation theory, one uses x = y? (3 / Attolr as a dimensionless 
expansion parameter, where an denotes the renormal- 
ized an {i-c. the denominator of the propagator). 
The zero-temperature limit corresponds to ax —* — oo or 
x — > oo. The parameter an is usually estimated via the 
Hartree approximation, where an — uh(1 — Lr) -1 M. 
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Using our non-perturbative method, we will be able to 
obtain a much improved estimate for x and olr. 

In the renormalization group analysis of the general- 
ization of (]]]) in 6 — e dimensions, Brezin et al. in- 
dicated that the GL hamiltonian in (|l|) is not closed 
under renormalization. The renormalization drives the 
quartic vertex of ([!]) to the most general gauge invari- 
ant form, /n j= i,2 d <^ e "^ (l " l|2+|Z2|2) l^( z i)| 2 5(Ni " 
Z2|)|0(22)| 2 , for an arbitrary function g. In this Let- 
ter, we use the Fourier transform representation || of 
the quartic vertex instead of the one used in Ref. ||: 
ffW — / d 2 r g(r) exp(i\i ■ r). It is convenient for 
later use to introduce the scaled function f(k) = 
(2//3) exp(— k 2 /2fi 2 )g(k). Note that the bare interaction 
corresponds to f R {k) — exp(— fc 2 /2^ 2 ). 

In the present work, we obtain the renormalized quar- 
tic vertex function, f R (k) of the GL hamiltonian by 
summing all the parquet graphs ||. We then calculate 
the structure factor of the two-dimensional vortex liq- 
uid from this non-perturbative solution. The structure 
factor denoted by A(q) is proportional to the Fourier 
transform of the spatial correlation function of the su- 
pcrfluid density, x(q) = / d 2 R e jqR x(r, r + R), where 
X (r,r')^<|*(r)| 2 |vl/(r')| 2 )-<|^(r)| 2 )(|*(r')| 2 ): 
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The structure factor of the two-dimensional vortex liq- 
uid has been studied using the numerical simulations [B, 
the density-functional theory of vortex liquid freezing pQ ] 
and the high-temperature perturbation theory ]TT[ ] . In 
terms of f R (k), A(q) is given by 



A(Q) = 1 - 2x{f R (Q) 



dKKJ (QK)f R {K)}, (3) 



where Jo is the Bessel function and Q = q//i and 
K = k//i are dimensionless wave vectors. Eq. (||) is the 
central quantity of our investigation. In (^), the first 
term comes from contracting two pairs of external lines 
in the disconnected four-point Feynman diagram. The 
remaining terms represent the two independent ways of 
contracting the connected four-point diagrams. 

Using f R (K), one can also express a R in terms of an- 
The relation between the renormalized parameter a R (or 
x) and an (or ax) gives the specific heat of the sys- 
tem. In the present analysis, this relation comes from the 
Dyson-Schwinger equation associated with the hamilto- 
nian ([!]), which expresses the renormalized propagator in 
terms of the bare propagator and the renormalized four- 
point function. It follows that 



1 - 4x{l -2x1 dK K f R (K)e' K /2 } 



The Hartree approximation therefore corresponds to 
neglecting the term that depends on f R (K). The 



generalized Abrikosov ratio is defined by (3a(x) = 
(|*(r)| 4 )/[(|*(r)| 2 )] 2 , where the bar denotes the spatial 
average. In terms of f R (K), we have 



Pa(x) = 2 



l-2x 
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It is interesting to see if /3a (x) really converges to 1.16 
as x — * oo or ax ^ — oo. 

In order to calculate f R (K), we consider the Feyn- 
man diagrams contributing to the connected four-point 
function, W(z\ , z|, z 3 , z 4 ) = (0*(^)0*(z|)0(2 3 )^ 4 )) c . 
We then make the parquet graph analysis on W and ob- 
tain the corresponding relations for f R (K) to which W 
is directly related within the LLL approximation 
The four-point vertex function depends on three variables 
characterizing the three independent channels. Since 
f R {K) is defined through the faithful representation of 
diagrams, we express all the diagrams involved in the 
parquet decomposition for W in the faithful representa- 
tion. 

Since an arbitrary diagram contributing W is either 
reducible or irreducible for a given channel i — 1,2 or 3, 
one can write in terms of f R (K) that f R (K) = h(K) + 
Ti(K), where Ii and r, are the irreducible and reducible 
parts for channel i, respectively. The irreducible part 
Ii(K) consists of the bare vertex function f R (K), the 
totally irreducible part R(K), which is irreducible in all 
channels, and those which are reducible in the other two 
channels. Thus, we have 

h{K) = f B {K)+R{K) + Y d Tj{K)- (5) 
Therefore, 
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f R (K) = f B (K) + R(K) + (6) 




FIG. 1. Diagrammatic representation of the Bethe-Sal- 
peter equation for the reducible parts Ti(K). 



2 



The reducible part Tf(K) can be represented in terms 
of the Bethe-Salpeter equation as described graphically 
in Fig. [I]. Essentially, they express the reducible part in 
terms of a series of ladder graphs composed of appropri- 
ate irreducible parts. The detailed derivation of these 
relations will be presented in a separate publication Jllj . 
By explicitly evaluating each diagram in Fig. [I], we obtain 

V 1 (K) = -2x J |^(|K-P|) (7a) 

x/ fi (F)cos(KxP) 
T 2 (K) = -2x{h{K)f R (K) (7b) 

+ h(K) J ^-fn(P) cos(KxP) 

+ f R {K) J ^/ 3 (P)cos(KxP)} 

r 3 (*Q = -2xj^ Z 3 (|K - P\)f R (P), (7c) 

where K x P = K1P2 — K2P1. It is important to note 
that the coupled integral equations, (||) with (||) and (0) 
are exact relations for f^(K). In the present work, we 
neglect the contribution from the nonparquet diagrams 
represented by R(K) in (|J), for which a systematic anal- 
ysis is very difficult. We note that the lowest order di- 
agram contributing to R{K) is of 0(/3 4 ). The subset of 
parquet graphs contain an enormous number of diagrams 
and we assume that the essential properties of the vor- 
tex liquid can be represented by this subset. We solve 
numerically eqs. (||), (||) and (Q) for fn{K) with the con- 
dition R(K) = and evaluate A (A) using (||). For given 
set of functions, {Ii(K)}, Eq. (||) is a linear equation for 
fii(K) (a Fredholm equation of the second kind) which 
can be solved by a numerical inversion of a matrix p3| . 
The solution is then found by iteration with the improved 
irreducible parts calculated from (||) using the solution to 
the linear equation just obtained. 

We obtain fn(K) and A(A) for various values of x. 
Before presenting our results for A (A), we first discuss 
how the perturbative expansion parameter x depends on 
dimensionless temperature ax- From our numerical data, 
we find a significant improvement over the Hartree ap- 
proximation. We also find that at low temperatures the 
yjx vs. ctT curve approaches a straight line, the slope of 
which is related to the Abrikosov ratio in the zero tem- 
perature limit. We can in fact plot directly the Abrikosov 
ratio as a function of temperature using (Q). (See Fig. ||.) 
The value of Pa decreases with decreasing temperature 
from its high temperature limit Pa = 2 representing the 
uncorrelated vortex liquid. As one goes into the low 
temperature region, Pa decreases very slowly. Although 
it is hard to determine the exact asymptotic value at 
qt = — 00 from the present data, these values are get- 
ting close to 1.16, the value for a triangular lattice. 
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FIG. 2. The generalized Abrikosov ratio as a function of 
temperature. 

Fig. H shows the structure factor A (A") for various 
temperatures. We find that, as olt —2, the struc- 
ture factor starts to develop a well-defined peak near the 
first reciprocal lattice vector (RLV) of a triangular lat- 
tice, Gi/fi ~ 2.69. As the temperature is lowered fur- 
ther, more structures are revealed as the peaks around 
the next RLVs begin to appear. Since the second and 
third RLVs of the triangular lattice are closely spaced, 
our solution could not resolve these peaks up to the min- 
imum temperature we have considered. However, one 
can clearly see that the height of the first peak grows 
while the width narrows down as the temperature is low- 
ered, which indicates the growing crystalline order in the 
two-dimensional vortex liquid system. 

In Fig. [5], our nonperturbative results for A (A) are 
compared with the recent MC simulation data H|. Up 
to low temperatures, they are in good agreement with 
the MC results. But we find that, as one goes further 
into the low temperature regime, the peak height in the 
MC results grows faster than in the present calculation. 
It is to be expected that the diagrams omitted from the 
parquet approximation will produce quantitative errors 
at low temperatures, whereas at high temperatures their 
effects are negligible. In Ref. pl| , the structure factor 
of the two-dimensional vortex liquid has been evaluated 
up to 12th order of perturbation theory. The extrapola- 
tion to low temperatures via Pade approximants shows a 
peak starting to appear near the first RLV of the trian- 
gular lattice, but the method becomes unreliable at low 
temperatures and at small wave vectors. 

In the completely ordered state where vortices form a 
lattice, the structure factor would exhibit delta-function 
like divergences at the reciprocal lattice vectors (RLV). 
In fact, if we use the mean field solution for ^(r) Q| 
in the definition of the structure factor, we get A(k) = 
(27r) 2 E G #o (5(2) (k-G), where the sum is over the RLVs 
of the triangular lattice except G = 0. The G = term 
corresponds to the disconnected piece substracted in the 
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FIG. 3. The structure factor A(K) calculated from the 
present non-perturbative method for various temperatures; 
olt = —1.5 (the broken line), —3.9 (the dotted line) and —10 
(the solid line) . The arrows indicate the positions of the RLVs 
of the triangular lattice. The open circles and the open trian- 
gles are those obtained in the MC simulation for qt = — 1.5 
and —4.0, respectively. The inset shows the height of the first 
peak obtained from the value of A at the first RLV, G\ as a 
function of temperature. 

definition of x- Assuming that A(K) obtained from our 
calculation approach this low temperature limit, we in- 
vestigate the temperature dependence of the height of the 
growing peaks in the structure factor. From our numeri- 
cal solution, we find that in the low temperature regime 
the height of the first peak grows as \ar\- (See the in- 
set in Fig. ||.) This result can be understood in terms of 
the zero-temperature scaling argument Q mentioned be- 
fore. This argument is based on the observation that the 
destruction of off-diagonal long range order below four 
dimensions |l5f| by thermal fluctuations might be accom- 
panied by the destruction of the positional order charac- 
terized by a length scale £. In this argument, £ diverges 
in the zero-temperature limit as \cxt\- If w ^ can describe 
A(fc) near a peak G by the function £h(£(k — G)), then 
the integral of A(fc) over k around G is finite as £ — > oo 
jjj]], which is required by the asymptotic form of A(fc). 
Therefore, the peak height at k = G grows as £ ~ \ar\ 
as otr — > — oo. 

We note that our non-perturbative approximation cap- 
tures this asymptotic behavior suggested by the zero tem- 
perature scaling. From (^), we deduce that in the zero 
temperature limit, /fl(k) ~ (2x) _1 {l — nfJ? J2a <5( 2 )(k — 
G)}. If we put this expression into eqs. (g), (||) and 
(fjj) and focus on the region k ~ G ^ 0, we find 
(4n/x)S^ (k — G) = a smooth function of k around G 
independent of x. Upon integrating over the region with 



temperature scaling and therefore within our approxi- 
mation the phase transition to the vortex lattice state at 
finite temperature is highly unlikely. 

To summarize, we have found an analytic approach to 
the vortex liquid regime which is sophisticated enough to 
predict, ab initio, the growth of crystalline order as the 
temperature is lowered; which agrees qualitatively with 
MC estimates of the structure factor at moderate to high 
temperatures; while in the low-temperature regime it is 
consistent with zero-temperature scaling argument, and 
the absence of a finite temperature phase transition. 

We thank T.Newman for useful discussions and 
M.Dodgson for providing his MC simulation data. 
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This argument indicates that our data for A(K) already 
exhibit the asymptotic behavior suggested by the zero- 
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